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For certain non-magnetic solids with low symmetry the occurrence of spin-polarized longitudinal 
currents is predicted. These arise due to an interplay of spin-orbit interaction and the particular 
crystal symmetry. This result is derived using a group-theoretical scheme that allows investigating 
the symmetry properties of any linear response tensor relevant to the field of spintronics. For the 
spin conductivity tensor it is shown that only the magnetic Laue group has to be considered in this 
context. Within the introduced general scheme also the spin Hall- and additional related transverse 
effects emerge without making reference to the two-current model. Numerical studies confirm these 
findings and demonstrate for (Aui-ajPtx^Sc that the longitudinal spin conductivity may be in the 
same order of magnitude as the conventional transverse one. The presented formalism only relies 
on the magnetic space group and therefore is universally applicable to any type of magnetic order. 


The discovery of the spin Hall effect [1-3] (SHE) with 
its particular feature of converting a longitudinal charge 
current into a transverse spin current has sparked nu¬ 
merous studies that finally led to a deep understand¬ 
ing of many effects that are spin-orbit induced. Among 
them are the enigmatic anomalous Hall effect (AHE) that 
shares the same origin as the SHE and many new phe¬ 
nomena emerging from a coupling of spin-, charge- and 
orbital degrees of freedom in electric fields as well as tem¬ 
perature gradients. Examples for these are the Edelstein 
effect (EE [4, 5]), as well as the spincaloritronic pendants 
to the SHE and AHE, namely the spin- and anomal¬ 
ous Nernst effects (SNE [6, 7], ANE [8, 9]), respectively. 
Many models have been formulated that aim to capture 
particular contributions to theses effects. For instance, 
the concept of the semiclassical Berry phase that can be 
determined on the basis of the band structure of perfect 
crystalline systems is connected to so called intrinsic con¬ 
tributions [10-12], Extrinsic contributions arising from 
scattering at impurities in non-perfect systems can, for 
example, be obtained from diagrammatic methods [13] or 
Boltzmann transport theory [14]. 

The aforementioned transport phenomena and their 
different contributions being linear in the driving fields 
should, in principle, be described using the funda¬ 
mental Kubo formula for the corresponding response 
function [15], 


Tij(u) = / d te~ lu,t / dATr pAjBi{t + \h\) 
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The effects then emerge from the characteristics of the 
underlying Hamiltonian, the pair of chosen operators for 
perturbation [Aj) and observable ( Bi ) and the symmetry 
of the system. Due to the intractability of the problem 
to exactly solve the Kubo formula for a realistic system 
in practice one has to recourse to approximations and/or 
models. However, irrespective of this problem one can 
still analyze the transformation properties of response 
tensors r determined by the Kubo formula to make state¬ 
ments which effects are in principle allowed, i.e. which 


non-vanishing tensor elements may occur given a particu¬ 
lar transformation property of the operators appearing in 
Eq. (1). This route has been followed by Kleiner [15, 16], 
who demonstrated that the occurrence of the AHE is pre¬ 
dicted by such a space-time symmetry analysis. Further¬ 
more, considering in addition heat currents he derived 
general Onsager reciprocity relations. 

Here, by extending this approach and applying it in 
the context of spin current operators [17] we demonstrate 
that in certain non-magnetic low symmetry systems an 
electric field can induce a longitudinal spin polarized cur¬ 
rent [18] that has hitherto evaded perception, and com¬ 
plements the transverse spin Hall effect. Furthermore two 
additional transverse effects are found which differ from 
the SHE by the direction of polarization. The results 
of the group theoretical analysis are independently veri¬ 
fied for an alloy bulk system performing relativistic first- 
principles Kubo-type transport calculations. The presen¬ 
ted formalism is furthermore very general, because (i) it 
allows to identify other response phenomena as non-zero 
elements in respective response tensors, as e.g. the AHE, 
(ii) it applies to both, magnetic and non-magnetic sys¬ 
tems, and (iii) it is free of the notion of a two-current 
model often used as an approximation in discussing spin- 
tronic phenomena; instead it is based on the concept of 
spin (polarization) current densities. 

The material-specific features of any transport prop¬ 
erty may be discussed on the basis of the corresponding 
response function tensor r. Concerning this, the shape 
of the tensor t, i.e. the occurrence and degeneracy of 
non-zero elements, reflecting the symmetry of the invest¬ 
igated solid, is obviously of central importance. To find, 
in particular, the shape of the spin conductivity tensor, 
Kleiner’s scheme [15] to deal with the symmetry proper¬ 
ties of ordinary transport tensors has been extended to 
the case when the response observable is represented by 
an arbitrary operator product of the form ( BiCj ) while 
an operator Ak represents the perturbation and the oper¬ 
ators Ak, Bj , and Ci are seen as the Cartesian compon¬ 
ents of vector operators. Within Kubo’s linear response 
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formalism the corresponding frequency (cu) dependent re- sponse function r is then given by 

I 


f°o w3 

T (B i c j )A k ( u ' H) = J Q dATr 


where as usual [15] p stands for the density operator, 
/3 = 1 /ksT with ks the Boltzmann constant, T is the 
temperature and H is a magnetic field that might be 
present. 

The shape of r can be found by considering the impact 
of a symmetry operation of the space group of the solid on 
Eq. (2), as this will lead to equations connecting elements 
of r. Collecting the restrictions imposed by all symmetry 
operations the shape of r is obtained. In this context 
it is important to note that the relevant space group of 
the considered system may contain not only unitary pure 
spatial ( u) but also anti-unitary symmetry operations (a) 
that involve time reversal. 

The transformation properties of the operators X = 
Ai, Bi or Ci in Eq. (2) under symmetry operations can 
be expressed in terms of the corresponding Wigner D- 
matrices [15] D^(u) and D {X \a) belonging to the op¬ 
erator X and the operation u or a, respectively. Starting 
from Eq. (2) and making use of these transformation re¬ 
lations one gets the transformation behavior of r under 
a unitary ( u ) or anti-unitary (a) symmetry operation, 
respectively: 

Imn 
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It should be noted that in general the tensors T^g.c-)A k 
and T 4 t (Btcf) are different objects representing differ¬ 
ent response functions which are only interrelated by 
Eq. (4) . It nevertheless imposes restrictions on the shape 
of giving rise to (generalized) Onsager rela¬ 

tions. 

Assuming Ci = 1 and Bi = Ai = ji with ji the current 
density operator r corresponds to the ordinary electrical 
conductivity tensor <r. Using the behavior of ji under 
symmetry operations [15], it turns out that only the mag¬ 
netic Laue group of the system has to be considered, that 
is generated by adding the (space) inversion operation I 
to the crystallographic magnetic point group [19]. The 
resulting shape of the conductivity tensor g_ is given in 
Table I for four different magnetic Laue groups. 


( / 9(H)i fc B J (t + mA;H)C' J (t + i? l A;H)) , (2) 


When considering the spin conductivity tensor its ele¬ 
ments tJij give the current density along direction i for 
the spin polarization with respect to the fc-axis induced 
by an electrical field along the j-axis. In this case the 
perturbing electric field is still represented by A, = ji 
while the induced spin current density is represented by 
the corresponding operator jf = ( BiC k )• As the expli¬ 
cit definition of Jf is not relevant for the following, but 
only its symmetry properties, the frequently used non- 
relativistic definition = ^{vi,(T k } may be used that 
consists in a combination of the Pauli spin matrix a k and 
the conventional velocity operator Vi [20]. Alternatively, 
one may use the relativistic definition of the spin current 
operator jf = Tkji as suggested by Vernes et al. [21] 
that involves the spatial part Tk of the spin polarization 
operator [22]. 

Expressing the transformation behavior of jA in terms 
of the Wigner matrices allows deducing the shape of the 
corresponding spin conductivity tensor on the basis of 
Eqs. (3) and (4). As for the electrical conductivity it 
turns out again that one has to consider only the mag¬ 
netic Laue group, i.e. there are only 37 different cases. 
Table I gives for the four cases considered here the shape 
of the various sub-tensors <r fe , where k specifies the com¬ 
ponent of the spin polarization. 

Considering a non-magnetic metal with fee- or bcc- 
structure (m3ml / ) Kleiner’s scheme naturally leads to an 
isotropic electrical conductivity tensor <x. The extension 
to deal with the spin conductivity tensor sketched above 
gives in this case only few non-vanishing elements that 
are associated with the SHE and are symmetry related 
according to: <r* z = crj x = = -<r* y = -a\ 7i = -cr^ x , 

i.e. cyclic permutation of the indices gives no change 
while anticyclic. permutation changes the sign. In con¬ 
trast to other derivations, there is obviously no need to 
artificially introduce a spin-projected conductivity nor 
to make reference to the conductivity tensor of a spin- 
polarized solid. For a ferromagnetic metal with fee- or 
bcc-structure (4 /mm'm') with the magnetization along 
the z-direction, the well-known shape of the conductiv¬ 
ity tensor <r is obtained that reflects the anomalous Hall 
effect (cr X y) as well as the magneto-resistance anisotropy 
(u xx 7^ o’zz) with the symmetry relations a xy = — cr yx and 
<7 XX = Uyy. As one notes the spin conductivity tensors 
show as for the non-magnetic case off-diagonal elements 
that represent the transverse spin conductivity. This im¬ 
plies the occurrence of the spin Hall effect in ferromag- 
nets that was investigated recently for diluted alloys [23] . 
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Table I. Electrical (a) and spin ( a _ k ) conductivity tensor forms for the magnetic Laue groups discussed in the text [18]. Below 
each group symbol an example for a material is given in parentheses. 


However, more elements show up as compared to the non¬ 
magnetic case since less symmetry relations survive in the 
presence of a spontaneous magnetization. Additionally, 
in contrast to the non-magnetic case also a longitudinal 
spin-polarized conductivity (cr^) occurs in a ferromag- 
net, that for example gives rise to the spin-dependent 
Seebeck effect [24]. A simple explanation for the corres¬ 
ponding longitudinal spin transport would be based on 
Mott’s two-current model assuming different conductiv¬ 
ities for the two spin channels. However, it is well known 
that spin-orbit interaction leads to a hybridization of the 
spin channels and influences even the longitudinal con¬ 
ductivity of a ferromagnet this way [25] . Accordingly, it 
cannot be ruled out that the longitudinal tensor elements 
crfj are not only reflecting the spontaneous spin magnet¬ 
ization of the material but are to some extent due to 
spin-orbit coupling. 

Indeed the scheme presented above leads for non¬ 
magnetic systems having low symmetry not only to off- 
diagonal elements reflecting transverse spin conductivity, 
i.e. the SHE, but also to diagonal elements reflecting lon¬ 
gitudinal spin transport, that was not observed so far. 
For the two magnetic Laue groups A/mV and 2/mV for 
non-magnetic solids considered in Table I, a 4- and 2-fold, 
resp., rotation axis is present. As a consequence longit¬ 
udinal spin currents show up only with spin polarization 
along this principal axis of rotation. 

To verify the results of our group theoretical approach 
independently we calculated the full spin conductivity 
tensor for solids having different structures corresponding 
to different magnetic Laue groups. This work employs a 
computational scheme that has been used before for nu¬ 
merical studies on the SHE in non-magnetic transition 
metal alloys [26] . Performing these calculations without 
making use of symmetry led numerically to a spin con¬ 


ductivity tensor that was always fully in line with the 
analytical group-theoretical results concerning the shape 
and degeneracies of the tensor. 

To get a first estimate of the order of magnitude 
of the longitudinal spin polarized conductivity in non¬ 
magnets, calculations have been done for the system 
(Aui_ x Pt x ) 4 Sc having the magnetic Laue group A/ml' 
for varying Pt concentration x. Fig. 1 (top) shows the 
corresponding electrical conductivity that is, in agree¬ 
ment with Table I, diagonal and slightly anisotropic, i.e. 
cr xx = cryy ~ <7 ZZ . Furthermore, the conductivities an are 
strongly asymmetric w.r.t. to the concentration x when 
replacing Au with prominent sp-character at the Fermi 
level by Pt with dominant d-character. Furthermore, one 
notes a relatively strong impact of the vertex corrections 
on the Au-rich side of the system (x ss 0) while these are 
much less important on the Au-poor side (x « 1). This 
observation is well known from binary transition metal 
alloys, like Cui_a;Pt x [27] or Agi-^Pdx [28], where the 
dominance of sp-character changes to d-character when 
x is varied from 0 to 1. 

The transverse spin conductivity <7* is shown in the 
middle panel of Fig. 1 for x-polarization of the spin. As 
Table I shows going from m3ml' to A/mV symmetry the 
relation <jy Z = — a zy disappears, i.e. the corresponding 
sub-tensor is not anti-symmetric anymore. A symmet¬ 
ric component, which is by definition not present in the 
ordinary SHE, indeed can be seen in Fig. 1 (middle) al¬ 
though the deviations are not very pronounced. In line 
one finds (except for x —> 0) for the additional non-zero 
tensor elements a* z « — a zx . The first coefficient relates 
a spin current jf c polarized in the direction of motion to 
an electric field E z , whereas <x zx describes a spin current 
j/ transverse, but with the spin polarization parallel to 
the driving electric field E x . To our knowledge the cor- 
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Figure 1. Top: longitudinal conductivity an for 

(Aui-rcPt x ) 4 Sc as a function of the concentration x calculated 
without (NC) and with (VC) the vertex corrections. Middle: 
transverse spin conductivities afj. Bottom: transverse and 
longitudinal spin conductivity a xy and <r xx , respectively. 


responding effects have not been considered so far. Inter¬ 
estingly, both elements occur simultaneously for a given 
magnetic Laue group or both are absent. However, com¬ 
pared to the spin Hall-like elements a yz and a * y they are 
smaller. For y-polarization of the spin the correspond¬ 
ing tensor elements are uniquely related to those for x- 
polarization according to Table I and for this reason not 
given here. The tensor elements afj for z-polarization are 
given in the lower panel of Fig. 1. In line with Table I they 
obey the symmetry relation a xy = —a yx (i.e. describing 
the pure SHE) and differ from a yz . This difference how¬ 
ever is, except again for x —>• 0 , not very pronounced. 
In particular, a yz and a xy show a similar variation with 
concentration x that differs clearly from that of the lon¬ 
gitudinal spin conductivity a xx shown as well in Fig. 1 


(bottom). Although this new type of tensor clement is 
overall somewhat smaller in magnitude than the dom¬ 
inating transverse elements it has nevertheless the same 
oder of magnitude, especially in the Au-rich regime, and 
for that reason it should be possible to determine it ex¬ 
perimentally. 

As can be seen in Fig. 1 the curves for the spin con¬ 
ductivity tensor elements cA as function of the concen¬ 
tration x are much more structured than the electrical 
conductivity an, i.e. they are much stronger affected by 
the variation of the electronic structure with composi¬ 
tion. In particular the spin conductivities afj show pro¬ 
nounced peaks or dips for x ~ 0.8. This behavior can 
be related to the variation of the density of states (DOS) 
with x as can be seen from Fig. 2. It shows the com¬ 
ponent resolved DOS n a (E ) as a function of the energy 
E for (Auo. 5 Pto. 5 ) 4 Sc (top) and at the Fermi energy Ep 
for (Aui_ x Pt x ) 4 Sc as a function of the concentration x 
(bottom). As mentioned above, at the Fermi energy the 



Figure 2. Top: energy dependent component (a) re¬ 
solved DOS n a (E) for (Auo.sPto. 5 ) 480 . Bottom: com¬ 
ponent resolved DOS n a (EF) at the Fermi energy Ef for 
(Aui-^Ptx^Sc as a function of the concentration x. 

partial DOS n Aa (Ep) of Au is dominated by sp-states 
while that of Pt has dominant d-character. The pro¬ 
nounced dip of the Pt DOS n Pt (Ep’) at x ~ 0.8 is appar¬ 
ently responsible for the prominent features in the spin 
conductivity curves shown in Fig. 1 (middle and bottom 
panels). 

As mentioned before, for the longitudinal conductiv¬ 
ity an inclusion of the vertex corrections has primar¬ 
ily an impact at the Au-rich side of the system. The 
same behavior is found for the transverse (cA-) as well as 
the longitudinal (cr^) spin conductivity components. For 


























5 


the transverse spin Hall conductivity it could be demon¬ 
strated that the contribution connected with the vertex 
corrections corresponds to the so-called extrinsic contri¬ 
bution that is primarily caused by the skew scattering 
mechanism [23, 26]. The very similar dependence of c>A 
and (J 7 j,j on the vertex corrections suggests that this ap¬ 
plies also for the longitudinal spin conductivity. 

In summary, a group-theoretical scheme has been 
presented that allows determining the shape of response 
tensors relevant for the field of spintronics. Application 
to the spin conductivity tensor gave a sound and model- 
independent explanation for the occurrence of the trans¬ 
verse tensor elements responsible for the spin Hall effect 
and two additional, closely related effects. In addition it 
was found that for low symmetry longitudinal elements 
show up in addition even for non-magnetic solids that 
were not considered before. Independent numerical in¬ 
vestigations confirmed these results and demonstrated 
for (Aui_ x Pt x ) 4 Sc that the longitudinal spin conduct¬ 
ivity may be in the same order of magnitude as the 
transverse one. It should be noted that the discussion 
of the spin conductivity tensor was referring to the de¬ 
limit lo = 0. However, the tensor forms given in Table I 
also hold for finite frequencies, implying the occurrence 
of the ac-counterparts to the discussed effects. In addi¬ 
tion, the formalism is applicable to numerous other lin¬ 
ear response phenomena as e.g. the AHE, anisotropic 
magneto-resistance (AMR), the Edelstein effect [4, 5] , 
Gilbert damping [29], spin-orbit torques [30], etc.. Fur¬ 
thermore, using the fact that the operators for electrical 
and heat currents share the same transformation proper¬ 
ties the presented formalism can be applied to spincalor- 
itronic phenomena as well. 
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